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Abstract. 

The ages and masses of neutron stars (NSs) are two fundamental threads that make pulsars accessible to other sub- 
disciplines of astronomy and physics. A realistic and accurate determination of these two derived parameters play an 
important role in understanding of advanced stages of stellar evolution and the physics that govern relevant processes. Here 

1 summarize new constraints on the ages and masses of NSs with an evolutionary perspective. I show that the observed P-P 
demographics is more diverse than what is theoretically predicted for the standard evolutionary channel. In particular, standard 
recycling followed by dipole spin-down fails to reproduce the population of millisecond pulsars with higher magnetic fields 
(B> 4 X lO'^G) at rates deduced from observations. A proper inclusion of constraints arising from binary evolution and mass 
accretion offers a more realistic insight into the age distribution. By analytically implementing these constraints, I propose 
a "modified" spin-down age (T) for millisecond pulsars that gives estimates closer to the true age. Finally, I independently 
analyze the peak, skewness and cutoff values of the underlying mass distribution from a comprehensive list of radio pulsars 
for which secure mass measurements are available. The inferred mass distribution shows clear peaks at 1.35 M^, and 1.50 
Mq for NSs in double neutron star (DNS) and neutron star-white dwarf (NS-WD) systems respectively. I find a mass cutoff at 

2 M0 for NSs with WD companions, which establishes a firm lower bound for the maximum mass of NSs. 
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EVOLUTION 

MSPs are spun-up second generation neutron stars that 
have accreted mass and angular momentum from their 
evolved companion during the LMXB phase [1, 2]. For 
these pulsars, I parametrize the transitional phase after 
the cessation of active accretion and the subsequent evo- 
lution. ^'^ 



This, in-turn, tunes the initial Pq/Pq fraction at birth 
(Figure 1, orange line). The dominant mechanism for 
energy loss (e.g., « = 3 for pure dipole radiation) shapes 
the diagonal evolutionary path (Figure 1, green line). 
The distribution of millisecond radio pulsars along the 
green diagonal line will be implicitly coupled with their 
Galactic birth rates {R). 



Parametrizing millisecond pulsar evolution 

The equilibrium period distribution (D) of millisecond 
X-ray pulsars (MSXPs) at the end of the accretion phase 
will be tightly constrained by the Keplerian orbital period 
at their Alfven radius (Figure 1 , cyan line). The geometry 
combined with the accretion rate (M) experienced during 
the LMXB phase ultimately limits the amount of angular 
momentum that can be transfered onto the neutron star 
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^ For details see [3, 4, 5] 

^ http://www.kiziltan.org/research.html: Full resolution color figures 
and movies are available at this URL. 



Observed vs. projected P-P distribution 

Synthetic millisecond pulsar populations that are con- 
straint only by standard evolution are sampled from the 
parameter space {D,M,R). In order to test whether the 
observed P-P diversity can be reproduced, I integrate the 
consistent populations along with populations that show 
the 'most diverse' P-P variations (Figure 2). The results 
can be summarized as follows: 

• No physically motivated Pq distribution has been 
able to produce the whole millisecond radio pulsar 
population consistently. No viable millisecond X- 
ray pulsar period distribution could subsequently 
produce the observed relative ratios of young/old 
pulsars with relatively high B fields. 
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FIGURE 1. Parameters that shape the P-P demographics of 
milhsecond pulsars. 



FIGURE 2. Projected P-P distribution of millisecond pulsars 
compared to observations. 



Observational selection effects cannot account for 
the rates at which higher magnetic field millisec- 
ond pulsars (e.g., PSR B1937H-21) are observed. 
The fraction of observed young/old millisecond ra- 
dio pulsars with high B fields is higher than what the 
standard model predicts by several orders of magni- 
tude. 

Millisecond radio pulsars with characteristic ages 
Tc > lO'*' yr are born with Pq « P and small spin- 
down rates. Figure 2 shows that the apparent enigma 
of millisecond pulsars with ages older than the age 
of the Galaxy is mainly a manifestation of very low 
accretion rates experienced during the late stages of 
their LMXB phase. 

It is necessary to posit the existence of a separate 
class of progenitors, most likely with a different 
distribution of magnetic fields, accretion rates and 
equilibrium spin periods, presumably among the 
LMXBs that have not been revealed as millisecond 
X-ray pulsars. 

A millisecond X-ray pulsar period distribution that 
has sharp multimodal features coupled with non- 
standard energy loss mechanisms may be able to 
reconcile for the joint P-P distribution of millisec- 
ond pulsars. 



AGE 

The standard method for estimating ages of pulsars is 
to infer an age from the rate of spin-down. The char- 
acteristic age (Tc = P/2P) is only an approximation of 
the generic spin-down age and may give realistic esti- 
mates for normal pulsars for which the initial spin pe- 
riods are likely to be much shoiter than their observed 



peiiods (i.e., Pq <C P). This method can fail, however, for 
pulsars with very short periods. For millisecond pulsars, 
the details of the spin-up process pose additional con- 
straints on the period (P) and spin-down rates (P) which 
may consequently affect the age estimate. 



A modified spin-down age for millisecond 
pulsars 

To estimate ages of millisecond pulsars, we propose 
the following method in which constraints arising from 
the recycling process and mass shedding are imple- 
mented: 




t(B>B,) = 



t{B<B,) 



where B^ is the locus of points extending from the in- 
tersection of the diagonal spin-up line and the vertical 
mass shedding limit. P^h is the mass shedding limit which 
is expected to be between 0.85-1.32 ms depending on 
the equation of state of the NS matter [6, 7]. Although 
a inherently has numerous sources of uncertainty, the 
corresponding minimum post-accretion period that de- 
fines the spin-up line depends on the total accreted mass, 
but is insensitive to other parameters [see, 8]. There- 
fore, we parametrically adopt the re-normalized coeffi- 
cient a = 2.6lo 4 ^ lO^s*''^ and use it as a fiducial value. 

While the modified spin-down age t gives a better es- 
timate for the true age than T^, it is still a strict upper limit 
(as is T( for MSPs). The projected underlying age distri- 



are mere underestimates. 
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FIGURE 3. Expected true age distribution for millisecond 
pulsars. Color represents the true age T,. Downward aiTows 
(I) are neutron star spin frequencies measured in LMXBs 
which are used as progenitor seeds to reconstruct the synthetic 
population. 

bution for standard evolution is reproduced in Figure 3. 
The conclusions are summarized as follows: 

• Characteristic age estimates are strict upper limits 
for millisecond pulsars. 

• A significant fraction of millisecond pulsars are 
born with apparent older ages. 

• Approximately 30% of the population is expected to 
be born below the Hubble line (t^ > lO'^yr). These 
millisecond pulsars must have experienced very low 
accretion rates during the LMXB phase. 

• There are two sources of age corruption: (1) Age 
bias, which is due to secular acceleration. Older 
millisecond pulsars may mimic younger ages be- 
cause the measured spin-down rates (P) may ap- 
pear to spin-down faster due to this effect, signif- 
icantly in some cases. This bias can be corrected 
for millisecond pulsars with well determined dis- 
tances and transverse velocities. (2) Age contami- 
nation, which is primarily driven by low accretion 
rates experienced during the LMXB phase. A sig- 
nificant fraction of the millisecond pulsar popula- 
tion is expected to mimic older ages due to this 
contamination. While a statistical assessment of the 
cumulative contribution of the "age contamination" 
(for the population) is possible [4], there are no triv- 
ial methods which can be used to estimate the con- 
tribution to an individual source. 

• We estimate that ~30% of millisecond pulsars have 
ages that are overestimated by a factor of two or 
more. As a consequence, the NS-NS merger rates 
based on the characteristic ages of binary pulsars 



• While a braking index of n=3 is consistent with the 
millisecond pulsar P-P demographics, it is challeng- 
ing to disentangle the combined effects of a possible 
early phase during which energy was lost by more 
efficient processes such as gravitational wave radia- 
tion (n=5), and a subsequent secular evolution with 
torques less efficient than dipole radiation (n<3). 

MASS 

In recent years, the number of pulsars with secure mass 
measurements has increased to a level which allows us to 
probe the underlying neutron star (NS) mass distribution 
in detail. Here, I briefly review the range of expected 
masses and how they compare to the mass distribution 
inferred from radio pulsar observations. A ciitical in- 
depth review with details on the numerical method can 
be found in [5]. 



Birth, accreted and maximum mass 

Masses of NSs at birth are tuned by the intricate details 
of the astrophysical processes that drive core collapse 
and supernova explosions [9]. A careful analytical treat- 
ment of the uncertainties that affect the Chandrasekhar 
mass gives 1.08-1.57 Mp as a reasonable range for the 
birth mass [5]. The expected amount of mass required 
to spin up pulsars to millisecond periods {Aniacc = 0.10- 
0.20 Mo) can be derived by angular momentum argu- 
ments [for details see, 5]. For the upper mass limit, mod- 
em equations of states predict a range between 2.2-2.9 
M,7,. It is still unclear, however, whether very stiff equa- 
tions of states (EOSs) that stably sustain cores up to the 
GR limit 3.2 Mq) can exist. 



Estimating the underlying mass distribution 

The methods used for NS mass measurements have 
intrinsically different systematics, and therefore require 
a more careful treatment when assessing the underlying 
mass distribution. There is no compelling argument to 
have a single coherent (unimodal) mass distribution for 
NSs that, we know, have gone through dissimilar evo- 
lutionary histories and possibly different physical pro- 
cesses that lead to their production [10]. With an ever 
increasing number of pulsars with secure mass measure- 
ments, we have reached an adequate number by which 
we can now test whether the implied mass distributions 
are consistent with different sub-populations. 
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FIGURE 4. Neutron star mas.s distribution. The underlying 
mass distribution of NSs show bi-modality with peaks at 1.35 
Mo and 1.50 Mq for NSs in DNS and NS-WD systems respec- 
tively. The solid lines show the central 68% probability range 
[as in 5]. 



The inferred distribution is shown in Figure 4 and the 
ramifications can be summarized as follows: 

• Neutron star masses show peaks consistent with 
an initial mass function at 1.35±0.14 and a 
recycled peak at 1.50±0.25 Mq. 

' Considerable mass accretion {Aniacc ~ 0.15 M,7,) 
has occurred during the spin-up phase. 

• A secure lower limit to the maximum mass of a NS 
exists at 2 M,^;. This excludes the softer equations of 
states as viable configurations for matter at supra- 
nuclear densities. 

• The lack of truncation close to the maximum mass 
cutoff suggests that the maximum NS mass is set by 
evolutionary constraints rather than nuclear physics 
or general relativity, and the existence of rare super- 
massive NSs is possible. 
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{pi, a^; data) = 

n ... K-f')^ 

Y\[2n{G^^S})Y''\ (3) 

as the parametric likelihood and 

i^{^o\data) = 

j J N{.//Q;ii,(7^)P{n,(7-\data)dixda^ (4) 

for the posterior inference, where for a mass distribution 

OT; = + w,-,/ = 1, ...,n (5) 

m is the pulsar mass estimate, ^ is the (unobserved) 
pulsar mass with an associated w error. For independent 
and Wi, this yields 

/(^,w) = N{.^;ix,a^)N{w;0,S^) 

= N{m;ii,a^+S^) (6) 

where the associated errors are w; ^ N{Q,Sf). 

Unlike conventional statistical methods, with a 
Bayesian approach it is possible to separately infer 
peaks, shapes and cutoff values of the distribution with 
independent reliability measures. This provides unique 
leverage to probe parameters which separately trace 
independent astrophysical and evolutionary processes. 
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